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We study a full causal bulk viscous cosmological model with flat FRW symmetries and where 
the "constants" G, c and A vary. We take into account the possible efl'ects of a c— variable into 
the curvature tensor in order to outline the field equations. Using the Lie method we find the 
possible forms of the "constants" G and c that make integrable the field equations as well as the 
equation of state for the viscous parameter. It is found that G, c and A follow a power law solution 
verifying the relationship G/c^ = k. Once these possible forms have been obtained we calculate the 
thermodynamical quantities of the model in order to determine the possible values of the parameters 
that govern the quantities, finding that only a growing G and c are possible while A behaves as a 
negative decreasing function. 
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^ ■ I. INTRODUCTION 

' In a recent paper (see [l]) we study a cosmological model with flat FRW symmetries filled by a perfect fluid and 
^~~> , where the "constants" G,c and A were considerate as function on time t. By different reasons exposed in such paper, 
■ we took into account the possible effects of a c— variable into the curvature tensor in order to outline the field equations. 

Through the Lie group method we studied the possible forms of the functions G and c that make integrable the field 
• equations. In this way we were finding that G and c follow a power law solution verifying the relationship G/c^ ~ k. 
^ [ But unfortunately we were not able to determine if G and c are growing or decreasing functions on time t. 

In order to determine if G and c are growing or decreasing functions we try to take into account thermodynamical 
' considerations, that is to say, we hope that thermodynamical restrictions help us to determine the behaviour of such 
, functions. For this purpose, in this paper we consider a cosmological model with flat FRW symmetries filled by a full 
' causal bulk viscous fluid and where the constants G, c and A are considered as functions on time t. Once we have 
lO . outlined the field equations (taking into account the possible effects of a c— variable into the curvature tensor) we 
' rewrite them in order to obtain a second order differential equation in order to apply the standard Lie procedure. 
' The study of this ode through the Lie group method allows us to obtain the precise form of the functions G and c 
O^' that make integrable the field equations as well as the equation of state for the bulk viscous parameter ^. 

As we will sec in section 3 the field equations only admit scaling symmetries (note that we are working with flat 
bJ[). FRW symmetries) i.e. we are studying a self-similar model. This fact obligates that G and c follow a power law 
solution, c = cot^^, verifying the relationship G/c^ = k, i.e. G and c are functions on time t but in such a way that 
this relationship must be verified. We find another restriction under this symmetry. The bulk viscosity ^, must follow 



the law ^ = k-yp^/'^, i.e. we have found a concrete equation of state for the viscous parameter 7 = 1/2. All these 
results arc in agreement with our previous paper ([2]). 

Once we have found the possible forms of the functions G and c we calculate the energy density finding in a first 
approach that the only physical solution imply that G and c must be a growing functions on time t since Ki > 0. 

In section 4 we will calculate all the physical quantities in order to complete the solution for our model. In particular 
we are interested in calculating the entropy in order to obtain some restrictions for the physical parameters as Ki 
and to elucidate if G and c are growing or decreasing functions. In this way we find an exact solution for A which 
behaves as a negative decreasing function. 

We end summarizing all these results in the last section. 

II. THE MODEL 

Following Maartens [:>] , we consider a Friedmann- Robertson- Walker (FRW) Universe with a line element 

ds^ = c^{t)df - fit) [dx^ + dy^ + dz^) , (1) 



'Electronic address: abelcal@ciccp.es 



2 



filled with a bulk viscous cosmological fluid with the following energy-momentum tensor: 

Tt = {p + p + ii) u.u'' - (p + n) 5'y , (2) 

where p is the energy density, p the thermodynamic pressure, 11 is the bulk viscous pressure and Ui is the four velocity 
satisfying the condition = 1. 

The gravitational field equations with variable G, c and A are: 



1 ff-^TTGlt), 



R^k - 7^.9-.fci? = -^77^T^k + A(t)g,fc. (3) 



Applying the covariance divergence to the second member of equation (3) we get: 



div ( ^77" + SfA ] = 0, (4) 



'■■^ " I c G j ' 8^G ' 



that simplifies to: 



p + 3(p + p)i7 + 3Fn = -— -P--4P-, (6) 

SttG G c 

where i7 stands for the Hubble parameter {H — f / f). Therefore, our model (with FRW symmetries) is described by 
the following equations: 

2H-2-H + iH^ = -^{p + Tl)^Kc\ (7) 
c 

3F2 = ^p + Ac^ (8) 

p + ^[p + p + U)H=-^-p^+Ap-, (9) 
sttG G c 




(10) 

We would like to emphasize that deriving (8) and taking into account (7) it is obtained (9), that is to say, the 
conservation equation (9) could be deduced from the field equation as in the standard case where the "constants" G, c 
and A are true constants (see [1]). 

In order to close the system of equations (7-10) we have to give the equation of state for p and specify T, ^ and r. 
As usual, we assume the following phenomenological (ad hoc) laws [3]: 

p^ujp, ^ = k^p^, T^Dsp^ T^ip-^ = k^p^-\ (11) 

where < < 1, and fc-^ > 0, Z?^ > are dimensional constants, 7 > and 5 > {5 — -^j^ so that < (5 < 1/2 
for < < 1) are numerical constants. Eqs. (11) are standard in cosmological models whereas the equation for r is 
a simple procedure to ensure that the speed of viscous pulses does not exceed the speed of light. These are without 
sufhcient thermodynamical motivation, but, in absence of better alternatives, we use these equations and expect that 
they will at least provide an indication of the range of possibilities. For the temperature law, we take, T = Dsp^ , 
which is the simplest law guaranteeing positive heat capacity. 

The growth of the total commoving entropy E over a proper time interval {to, t) is given by [3]: 



S(i)-E(to)=-T^ / 

J to 



rfi, (12) 



T 



where /c^ is the Boltzmann's constant. 

Therefore, with all these assumptions and taking into account the conservation principle, i.e., div{T^) = 0, the 
resulting field equations are as follows: 
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2ij + 3i/2_2-i/ = -^(p + n)+Ac2, (13) 

3ff' = ^p + Ac^ 14 

p + 3(ij + l)pi7 = -3i^n, (15) 

n+^.-3,i/-ln(3i7-Ty0, (17) 



where W = I I is a numerical constant. 



III. LIE METHOD 



The field equations (13-17) represent a system of odes with 5 unknowns. In order to integrate them and to obtain 
a complete solution for the proposed model we will need to make some simplifying hypotheses (this is always a 
dangerous way) or to study if the system admits symmetries. In this way, studying the possible symmetries that 
admit the system, we will be able to determine the possible forms for which such system is integrable. 

In order to apply the standard Lie procedure (see for example [■5]-[7]) we need to rewrite eq. (13-17) by a single 
one and imposing that such ode (second order ode) admits a concrete symmetry we will be able to determine the 
equation of the state for the viscous parameter i.e. ^ = k-yp'^ and the exact of the "constants" G and c for which the 
ode is completely integrable. 

We start with the assumption 11 = xp, with x e R~ (see [2] for a complete discussion of this hypothesis). The bulk 
viscosity evolution equation can then be rewritten in the alternative form 

+ k-^p"-^ = -3f3H, (18) 
P ^ 

where /3 = (i + 1) and (5 = (l - -f ) . 

Taking the derivative with respect to the time of this equation and with the use of the equation 

H = ^F-47ra^p, (19) 

obtained from the field equations (13) (where a = {1 + ui + x)). and taking into account eq. (15) 

H = -^- (20) 
6a p 

we obtain the following second order differential equation describing the time variation of the density of the cosmo- 
logical fluid: 



p=^-Asp'p + B-p + D%p^, (21) 



p c 



where A = s (1 - 7), B = ^ and Z) = 

We go next to apply all the standard procedure of Lie group analysis to this equation (see [(>] for details and 
notation). In this way we find that the admissible symmetries for our ode being determinate by the following system 
of pdes 

P^PP + ^p = 0, (22) 

2p' Qr/p, - + (^Asp' - b'-^ - pr^p + r, = 0, (23) 

2vtp - c-^^tt - ^vtp-' + Ash^p'-' + (^sAp' - B^^ +^b(^-^ + '^^- Sp^D^^p = 0, (24) 

m + Vt (sAp^ - + D^p' (-26 + Vp^ -^p-'v) D^p'i ( g - 2^ ] = 0, (25) 
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Solving (22-25), wc find that 



witli the constraints, from eq. (24) 



and from eq. (25) 



G c c2 

being a and b numerical constants and we will assume that k > 0. 
In order to solve (27), wc consider the following cases. 

1. Case I. Taking a = 0, 5 ^ 0, we get 

^'2 



c at + 26 



G ^ c G ,^^\ 

= 2-, =^^ = K (28) 



c" = — , =^ c{t)^K2e''^\ (29) 
c 

where the {Kif^_^ are integration constants. Therefore the equation (21) yields: 

p=^-^^p + BKip + Dnp', (30) 
p z 



the solution obtained through invariants is: 



— = — p « const. (31) 



this solution seems not to be physical. 



2. Case II. Taking 6 = 0, a 7^ 0, we get that the infinitesimal X \s X = tdt — "Ipdp which is precisely the generator 
of the scaling symmetries. Therefore the solution will be the same than the obtained one with the dimensional 
method. 

With these values of a and h equation (27) yields 

^ c(t) = i^2t^S (32) 



C'2 c' 



C t ' 

as we expected, c{t) follows a power-law solution. Therefore equation (21) yields 



P=^-^^PP + B^P + Dkp\ (33) 



the solution obtained through invariants is: 

P « i , p = Pot , (34) 



where po is a numerical constants that must verifies eq. (33) in such a way that 



1 I A(A+ JA^ %Dk (1 4- BA'i) + 4i:iK (1 -\- BKi) 

^^-2[- I' 

now, since pq must be a positive numerical constant, in a first approach, and taking into account the definition 
of the constants A, B and D we see that po > iff (1 -I- BKi) < 0, which imply that A'l > since A > 0, B < 
0, K > and D < (since >c < 0). Therefore A'l must be a positive numerical constant. 
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Once again if one insists in solving equation (33) through DA (applying the Pi-theorem) it is found that with 
respect to the dimensional base 5B = {p, T} each quantity has the following dimensional equation [p] = p, [t] — T 
and [A\ = p^^T^^ . Therefore, we find in a trivial way that: 





P 


A 


t 


p 


1 


-1 





T 





-2 


1 



3. Case III. Taking a, 6 ^ 0, we get 



hence equation (21) yields: 



^'2 



ac 



P ■ 



at 



P 



1 



c{t) = K2 {at + by 



A 



2^^ + ^^^ + ^"^'' 



the solution obtained through invariants is: 

dt 
T 



dp 
V 



As we can see case II is a particular situation of case III. 



(36) 



(37) 



(38) 



(39) 



Through the symmetry analysis of eq. (21) we have been able of determining the equation of state that must follow 
the viscous parameter ^, obtaining that 7 = 1/2, as well as, that under the scaling symmetry c and G must follow a 
power law solution in such a way that they must verify the relationship G/c^ = k. 



IV. A COMPLETE SOLUTION 

Once we know the behaviour of the quantities p, G and c, in this section, we are going to determine the rest of the 
quantities i.e. the scale factor and mainly the behaviour of the entropy with the hope that this quantity help us to 
determine the possible values of the constant Ki. For this purpose we begin considering eq. (15) 

p + 3{uj + l + >t)pH = 0, (40) 

which trivially leads us to the well known relationship between the energy density p and the scale factor / 

P = A^r'("+'+") or p^A^r^^, (41) 

where a = {lu + 1 + m:) . Now, taking into account the eq. (17) and simplifying it, we obtain, 

thereby obtaining p = p(t). On simplifying further, we obtain, 

l_ _ 

where 



P = Pot-', (43) 



K=[l-^- — ~^], po = {2Kk,f. (44) 
' 2 Ka 2a' 



From equation (41) we obtain: 



A^ 



l/3a 



f=['-^f] , i.e. /cxf^r™^. (45) 
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An important observational quantity is the deceleration parameter Q = ^ {^) ^ ^- The sign of the deceleration 
parameter indicates whether the model inflates or not. The positive sign of q corresponds to "standard" decelerating 
models whereas the negative sign indicates inflation. In our model, the deceleration parameter behaves as: 

q = -l + Y- (46) 

In this way we find that for w = 1 , >ir = —4/3 is a critical value since q = 0, therefore and under these considerations, 
q<Oif >c< -4/3 and g > if x > -4/3. 

We proceed with the calculation of the other physical quantities under the condition 7 = 1/2: 

^ = k^p'' o, {pot-y ~ k?^t-\ (47) 

T = Dsp^ = Ds{pot-^Y ^t-^ with S=^^, (48) 

T = ^p-^ = k^ {pot-^y^^^ , i.e. T = {2Ky^ t (49) 

We see from r = {2Ky'^ t that this result is in agreement with the theoretical result obtained in [.)]. For viscous 
expansion to be non-thermalizing, we should have t < t, or otherwise the basic interaction rate for viscous effects 
should be sufficiently rapid to restore the equilibrium as the fluid expands. 
The commoving entropy is 



E(t) = Eo / x^^+^-^dx = (cj + 1) Eo<"- (50) 

a J 

where Eq = ^ — jj-j — - — and we find that x > — (cj + 1) i.e. x e (—2, 0] . Wc notice that the parameter x weakly 
perturbs the perfect fluid FRW Universe. When x ~ Q, the commoving entropy assumes a constant value and we 
recover the perfect fluid case. 

Finally, we will use equation (16) to obtain the behaviour of the cosmological "constant" A. Since we know the 
behaviour of the constants G, c and p i.e. 

G = KC^ c(t) = AV^S and p = pot-^, (51) 

then eq. (16) yields 

^-2^=0, (52) 

STTKp C 

and substituting eq. (51) into eq. (52) it is obtained the following ODE 

A^mnnpo^t-'"''' =^ A^-8^/^^-^t-2(^^+i) (53) 
K2 A2 Ai + 1 

with ATi > -1. 

We would like to emphasize that for ATi > then A < and for ATi < then A > but unfortunately we do not 
able to fix a better value for the constant Ki. Finally it is observed that 

A^^. (54) 

as it is expected in this self-similar solution. 

As we have seen, with the followed way, we have not been able to determine a better value for the constant Ki 
than the obtained one in eq. (35). for this reason we try a new tentative which consists in considering the effects of a 
c— variable into the scale factor in such a way that when we calculate the entropy such effect may by reflected better 
than in our previous way. 



A. A new tentative. 



In order to take into account the effects of a c— variable into the scale factor we will determine it from eq. (14) 
since we already know the behaviour of G, c, p and A. Therefore, after a simplification eq. (14) yields. 



3i/2 ^ Sn^po ( 1 - JJ^^, ] (55) 
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therefore 



f = V I ^-( 'T\/, S '\ (56) 



finding again that K\ > —1. 

In this way the entropy formula (12) yields: 



3 Ki + 1 



E(t) = Eo^^^+'^"^ (57) 



where Sq = — ><\/6po 2Tnipo+{K^^r){s^ ■ niake 27rKpo = 1 (a free hypothesis), S(i) « Sq^^^ ^^^^ "^^-^ find 
therefore A'l > —uj, but we are not able to find a better value (or limit) for A'l. 
Or taking into account the field eq. (13) 

2H + 3H^ - 2—H = -Sttkpq ( {lu + x) + J'^^ ] (58) 
t \ Ai + 1/ 

which is a Riccati ode and which particular solution (scaling solution) is: 

H^j f^Kff^ (59) 

where 



(A'l + 1)" ± J (A'l + 1) - 3A f a (A'l + ly + Kf + A'l 
a = ^- . (60) 

finding that A'l > — 1, with a = {lo + A ~ Sttkpq and imposing the condition. 

(A'l + 1)"^ - 3A (a (ATi + if + Kf + A'l) > (61) 
(note that a may take negative values) therefore the entropy behaves as: 



As we can see following these ways we have only been able to determine that A'l > — 1 or that 
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(ifl+l)2+^(Jfl + l)<l-3A(a(jri + l)2+A'2^Kl) 



(63) 



LO+l 

but as "(/"(ATi) depends on the factors like A and a it is very difficult to determine a good value for this constant. 

V. CONCLUSIONS. 

In this paper wc have studied a fiat FRW cosmological model filled with a bulk viscous fluid and where the 
"constants" G, c and A are considered as functions on time t. In order to outline the field equations we have taken 
into account the possible effects of a c— variable into the curvature tensor. In this way we have been able to deduce 
the general conservation principle from the field equations obtaining an equation that coincides with the obtained 

one from the divergence of the right hand of the field equation i.e. div (^§[T^ + SjA^ ~ 0. Using the Lie group tactic 

we have been able to determine the equation of state for the bulk viscous parameter, ^ = kyp'^ , as well as the exact 
form for the "constants" G and c which make integrable the field equations. 

We have obtained that our model is self-similar i.e. admits a scaling symmetry iff 7 = 1/2 and c follows a power law 
solution i.e. c = cot^'^'^, while G and c must verify the relationship G/c^ = k. Under these results we have calculated 
the behaviour of the energy density finding, in a first approach, that the physical solution is only possible if ATi > 0, 
that is to say, G and c are growing functions on time t. We have also calculated the rest of the main quantities but 
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under physical considerations we have only been able to determine that Ki > —1. The latter result does not contradict 
our previous result Ki > 0. 

Therefore wc conclude that G and c are growing functions while A is a negative decreasing function. 
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